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Medium-Energy Ion Scattering study of strained holmium silicide nanoislands grown 

on silicon (100) 
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(Dated: July 29, 2008) 

We have used medium-energy ion scattering (MEIS) to quantitatively analyse the structure of 
holmium silicide islands grown on the Si(lOO) surface. Structure fitting to the experimental data 
unambiguously shows that the tetragonal silicide phase is present and not the hexagonal phase 
which is associated with the growth of nanowires at submonolayer coverages. Islands formed with a 
lower holmium coverage of 3 ML are also shown to be tetragonal which suggests that the hexagonal 
structure is not a low coverage precursor to the growth of the tetragonal phase. MEIS simulations 
of large nanoislands which include the effects of lateral strain relief have been performed and these 
compare well with the experimental data. 
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Whilst there have been many studies of rare-earth (RE) 
silicides on the Si(lll) surface, relatively little was known 
about the growth mechanisms on the Si(lOO) surface un- 
til the discovery of self-assembled nanowires by Preines- 
berger et al. [l| These novel structures form when a 
suitable RE metal, eg. Gd,^ S,J1 [1 [SB Sc [iDy, 
[H, M i, M Sm, i Er, 0, a M, Elll Ho @,lll or 
Y, [i^ is deposited onto a clean Si(lOO) substrate held 
at an elevated temperature. Characteristic 'wires' mea- 
suring up to a micrometre in length, and typically only a 
few nanometres wide were observed by scanning tunnel- 
ing microscopy (STM). This discovery and the potential 
technological applications of such conducting nanowires 
has motivated considerable interest over the last ten years 

il 
23 



iiu ya [131, M \m M, yj, m \m 

|24|, |25|, |26|, [23, [2^ . Chen et al. have demon- 
strated that the growth of nanowires with extremely high 
aspect ratios is the result of the anisotropic lattice mis- 
match that results from the growth of the hexagonal, 
defect- AIB9 RE sihcide on the Si(lOO) surface (see Fig.[l|) 

Depositing more RE metal causes islands to form with 
the surface displaying a c(2 x 2) periodicity. The struc- 
ture of these RE silicide islands has been proposed to 

or orthorhombic 
All three of 

these phases are known to exist in the bulk and their lat- 
tice constants have been measured and it is known that 
the bulk orthorhombic phase is only a small distortion 
of the tetragonal phase [231 ■ ^ side-view schematic of 
the ThSi2 structure is shown in Fig. [21 Evidence for the 
growth of the tetragonal form has generally been inferred 
from lattice mismatch arguments, along with STM mea- 
surements of step heights in the silicides on the Si(lOO) 
surface [Tjl • SXRD measurements also yield information 
which appears to confirm that Er silicide islands adopt 
the tetragonal form, with the surface c(2 x 2) periodicity 
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FIG. 1: (a) A top down schematic illustrating the lat- 
tice matching of the hexagonal defect-AlB2 structure to the 
Si(lOO) surface. The anisotropy in lattice mismatch creates 
the high aspect ratio nanowires commonly observed for several 
RE metals, (b) A side view schematic of the AIB2 structure. 



interpreted as being due to Si adatoms [12] ■ 

STM experiments reveal that the island morphology 
is very sensitive to the RE metal deposited and to the 
annealing temperature used. Those RE metals (Nd. Sm, 
Yb) that have a low anisotropy in their lattice match to 
the substrate form compact 3D islands [Sj. Those that 
have a high anisotropy form both elongated and com- 
pact islands, depending upon the growth temperature. 
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FIG. 2: A side view schematic of the tetragonal ThSi2 struc- 
ture. The orthorhombic GdSi2 structure of HoSi2 is essen- 
tially the same but with a = 4.03 A and b = 3.94 A, a differ- 
ence of just 2.25 %. 
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FIG. 3: View from above of the bulk terminated Si(lOO) sur- 
face showing the two incident beam directions used in the 
MEIS experiments. For the [Il0]/[ll0] geometry the an- 
gle with respect to the sample normal was 45° and for the 
[lll]/[2ll] geometry the polar angle was 54.74°. The scat- 
tering plane through the tetragonal silicide is shown for the 
[110]/[ll0] geometry in figureH 



Growth at lower temperatures (600°C) causes the for- 
mation of elongated islands and growth at temperatures 
above 650°C results in compact 3D islands. However, the 
two island morphoplogies have been observed to coexist, 
especially when using intermediate annealing conditions. 

Dysprosium silicide provides an interesting case study. 
The elongated islands have typical dimensions 2 nm high, 
15 nm wide, and 500 nm long and the smaller compact is- 
lands have typical dimensions 5 nm high, 50 nm wide and 
200nm long [22|. Using high resolution cross-sectional 
TEM, Ye ei al. demonstrated that elongated and non- 
elongated islands coexist when Dy silicide is grown at 600 
°C - 650 °C on the Si(lOO) surface [13 . They proposed 
that the structure of the elongated islands was hexago- 
nal, and the excessive stress within this structure (due 
to the large c-axis mismatch) is relieved through dislo- 
cations and tilting across the width of the island. The 
non-elongated islands were found to be a tetragonal or 
orthorhombic structure with only a small amount of tilt- 
ing required to relieve the stress, since the lattice match 
is relatively small in both directions. It was also noted 
that there was some expansion of the c-axis to relieve the 
stress in much the same way as has been demonstrated 
for the 2D and 3D RE siheides on Si(lll) [33,[3l|. How- 
ever, the recent HR-TEM study by He et al. claims that 
the compact 3D Dy silicide islands are in fact a fully re- 
laxed, and hence stress free, hexagonal form, whilst the 
elongated nanowire islands are tetragonal/orthorhombic, 
with the faulted stacking relieving the stress in the struc- 
ture [22i |. The delicate energy balance within the sili- 
cide in this particular example is also highlighted by the 
fact that there is both a tetragonal and an orthorhombic 
phase of Dy silicide in the bulk. 

Despite the many studies that have been conducted. 



it has been noted in recent work that no ideal crys- 
tallographic determination has been possible to clearly 
prove once and for all the true nature of the 3D islands 
[23 . [26l . [28j created when the rare earth metals form a 
silicide in this way. MEIS is of particular value in the 
study of these nanostructures since the large mass con- 
trast of Ho and Si allows the elucidation of the silicide 
structure free from substrate effects. This also means 
that the technique is able to isolate the regions of inter- 
est on the surface by only selecting the Ho signal for the 
structural optimisation. 



I. EXPERIMENT 

All MEIS data were obtained at the UK MEIS facility 
at STFC Daresbury Laboratory. The experiments were 
conducted under ultra- high vacuum (UHV), with typical 
base pressures of 2 x 10~^° mbar. Clean Si(lOO) samples 
(cut from a lightly doped n-type wafer) were prepared 
by e-beam heating to ~ 1200 °C for 1 min, then slowly 
cooled (< 100 °C min-i) between 1000 °C and 600 °C 
to ensure an ordered surface was obtained. Temperature 
measurements were made using an infrared pyrometer. 
The heating cycles were repeated until a sharp 2x1 low- 
energy electron diffraction (LEED) pattern was observed. 
The 3D silicide islands were then formed by depositing 
approximately 6 ML of Ho from a quartz crystal cali- 
brated tantalum boat source onto this clean Si(lOO) sur- 
face, which was held at ^ 650 °C during deposition and 
for 5 min afterwards. When a c(2 x 2) LEED pattern 
was observed and the Auger electron spectrum showed 
the samples to be free from contamination, they were 
transferred under UHV into the ion scattering chamber. 

The ion scattering data were taken using 100 keV H+ 



ions incident upon the sample and the scattered ions were 
detected using an angle-resolving toroidal-sector electro- 
static ion-energy analyser and its microchannel plate de- 
tector. The MEIS spectra confirmed that the samples 
were free of contaminants and data were acquired with 
a total dose of ^ 10^^ ions cm^^. Further details about 
the Daresbury MEIS facility can be found in the liter- 
ature [30, [22, |3J, [SJ, [33 . During each experiment two 
different incident beam directions onto the sample were 
used. These beam directions are shown with respect to 
the substrate in fig. [21 The notation [m]/[oMi] for each 
geometry defines the ingoing crystal direction [in] and an 
outgoing crystal direction [out] that lies in the detected 
angular range of the scattered ions. 



II. STRUCTURE FITTING: TETRAGONAL OR 
HEXAGONAL? 

Due to the contention in the literature regarding 
whether these islands take the hexagonal, tetragonal or 
orthorhombic form, MEIS proves particularly powerful 
since it can clearly demonstrate the presence of one of 
these structures. In this work we have used the XVegas 
code [36| which uses Monte-Carlo methods to simulate 
the blocking curves of a proposed structural model. 



Fig. |4(a) shows the simulated blocking curves obtained 
when differing numbers of layers of RE are present in the 
hexagonal structure for the [110] /[1 10] geometry. The in- 
terlayer separation of all the Ho layers has been optimised 
to give the best visual agreement with experimental data, 
with the two intermediate Si layers lying 1/3 and 2/3 of 
this layer separation from the Ho layer beneath. A Ho- 
Ho layer separation of 3.28 A has been found to give the 
best-fit. However, the quality of the fits obtained is very 
poor, clearly demonstrating that this is not the correct 
structure. To make this clearer, the fit for 6 layers of Ho 
is shown in Fig. 4(b) At scattering angles of 88 , 90 , 
100 , 106 , 118 , 127 and 142 there arc major discrep- 
ancies between the experimental data and the best-fit 
simulations. The fully relaxed hexagonal structure that 
was proposed by He et al. [22| does not explain the exper- 
imental blocking curves either, as such a structure would 
produce almost identical blocking curves, but would yield 
a smaller c-axis of the silicide. 

The tetragonal structure has also been fitted to the 
data, since it is assumed that the silicide takes the same 
a and b values as the Si substrate to form an epitaxial 
overlayer. This goes against the available data in the lit- 
erature regarding the structure of the bulk Ho silicidcs, 
as no tetragonal phase has been observed (though the 
orthorhombic GdSi2 form has been reported with lattice 
parameters of a = 4.03 A, 6 = 3.944 A and c = 13.30 A) 
[291 ] . Conducting simulations of such a tetragonal struc- 
ture gives the fits for the [110]/[110] geometry shown in 
Fig. 5(a) which are for structures that are one, two and 
three ThSi2-type cells in height. It is clear that a crys- 
tal with a depth of two ThSi2 cells produces the best 



Experimental data 

2 Ho layers 

3 Ho layers 

4 Ho layers 

5 Ho layers 

6 Ho layers 






'\.t/\N 




- 2000 

- 1500 

- 1000 

- 500 



80 90 100 110 120 130 140 150 
Scattering angle (degrees) 
(a) 




80 90 100 110 120 130 140 150 
Scattering angle (degrees) 



(b) 



FIG. 4: Original in colour, (a) A comparison of the simula- 
tions of the hexagonal structure with 2-6 layers of Ho, to the 
experimental data for the [110]/[ll0] geometry, (b) The fit 
obtained for 6 Ho layers in the hexagonal structure, with a 
Ho-Ho layer separation of 3.28 A. 



match. The detailed structure fits are shown in figures 
[5(b)] and 



5(c) 



The r.m.s. thermal vibrations of the Ho 
atoms were set to the bulk metal value of 0.13 A, whilst 
a fitting of the Si vibrations gives an enhanced value of 
0.15 A. These vibrations are consistent with the previ- 
ous MEIS studies of 2D and 3D RE silicides on Si(lll) 
[31I . [SSj , with the factor of two enhancement for the Si 
atoms possibly indicating the presence of static disorder 
in these layers. This could be due to the presence of Si 
vacancies [20, 123] . 

Optimising the layer separation of the Ho atoms gives 
the best visual agreement for both geometries when c = 
13.14 A. It is clear from comparing fig. 5(b) with fig. [4(b)[ 
that the tetragonal form of the silicide is a better fit to 
the data than the hexagonal structure. 



III. AN ANOMALOUS C-AXIS VALUE 

The assumption that the silicide grows epitaxially on 
the Si(lOO) surface yields a constant c-axis value through- 
out this silicide (13.14 A) that is much smaller than the 
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FIG. 5: Original in colour, (a) Comparison of simulations of 
the tetragonal structure with thicknesses of 2-4 ThSi2 cells in 
height to the experimental data for the [110]/[110] geometry, 
(b) The fit obtained for the [IlO]/[110] geometry when the 
silicide is two ThSi2 cells in height, (c) The fit obtained for 
the [111]/[211] geometry when the silicide is two ThSi2 cells 
in height. 



value of 13.30 A for the bulk silicide. We would expect 
that the stress induced through a contraction of the Ho 
silicide a- and 6-axes at the interface, would be released 
through an expansion in the perpendicular direction and 
such behaviour has been directly observed in other rare 
earth silicide surfaces [30, l3l|. In particular, Ye et al. 



have reported an expansion in the c-axis of the silicide 
relative to the bulk value for Dy on Si(lOO) [131 . 

One explanation of this observation could be that the 
silicide is incommensurate with the Si substrate. Hence 
it would seem reasonable to assume that the silicide takes 
the orthorhombic form that is observed for bulk Ho sili- 
cide. Refitting the experimental data to simulations of 
this structure, where a = 4.030 A and b = 3.944 A, yields 
almost identical fits, but with a c-axis of 13.70 A. This 
large value of the silicide c-axis does not make sense ac- 
cording to the established trend in strain either, as it 
would be expected that the lack of strain at the interface 
would not cause any expansion in the c-axis of the silicide 
relative to the bulk value. Thus it can be concluded that 
it is not the bulk form of the silicide that is present on 
the surface; it must be some intermediate form. 

One of the limitations of MEIS when only using the Ho 
blocking curves to solve crystallographic structure is that 
it is not possible to exclusively determine both the a- and 
c-axis lattice constants of the silicide. This is because any 
difference in the silicide a-axis could be compensated by 
a change in the c-axis to produce blocking curves with 
blocking dips in the same positions. Hence, a particular 
c/ a ratio will yield a series of blocking curves that are 
exactly the same, even though the actual lattice parame- 
ters may be very different. We have performed a series of 
simulations, each with a fixed a-axis, in which the c-axis 
has been optimised by varying the layer separation of a 
silicide that is two ThSi2 cells in height. These results 
are shown in Table [J 
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TABLE I: Optimised c-axis values determined for various 
fixed values of the a-axis. Fitting the [ii0]/[ll0] and 
[lli]/[2ii] geometries yield slightly different values, which 
can be averaged to give an overall best-fit value. 



The data from Table |T] are plotted in Fig. [S] and the 
straight line fit defines a c/a ratio of 3.42 ± 0.01. The 
positive gradient suggests that as the lateral lattice con- 
stant of the tetragonal unit cell is reduced to match that 
of the substrate, the lattice constant perpendicular to the 
surface also contracts and this is true for all pairs of c, 
a values. As noted earlier, we might expect the opposite 
to occur and the c-axis to expand in order to maintain 
the volume of the unit cell. 
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FIG. 6: A plot showing how the c-axis determined is depen- 
dent on the size of the a-cixis of the silicide. Fitting a straight 
line which passes through the origin yields a c/a ratio of 3.42 
± 0.01. 



IV. SIMULATION OF STRAINED ISLANDS 
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FIG. 7: View from above of the 10x10x2 unit cell showing 
the lateral strain. 



There is another reason that might explain the unusu- 
ally low c-axis suggested by the fitting procedure. If the 
lateral lattice constant of the unit cell is constrained to 
match that of the substrate at the interface and is al- 
lowed to expand to relieve strain as we move away from 
the interface then the blocking angles will be shifted. The 
average a-axis throughout the structure will increase and 
so will the fitted c-axis in order to maintain the c/a ratio 
of 3.42. This effect would allow for interfacial matching 
and would at the same time provide a more physically 
meaningful thickness for the silicide. 

Simulation of this effect is computationally intensive. 
Unit cells with lateral strain as a function of height do 
not laterally tesselate ad-infinitum and a structure with 
finite lateral extent must be considered. This must be 
large enough to approach the size of a typical surface is- 
land but small enough to be computationally tractable. 
We have taken a unit cell with a 10x10 lateral extent, 
making an island 38.4x38.4 A^. This is two tetragonal 
unit cells deep and in total includes 2500 atoms. Com- 
putational resources prohibit us from a full optimisation 
of the strained structure in terms of interlayer spacings 
or vibrational amplitudes. Instead we have taken val- 
ues suggested from the structure fitting carried out on 
the 1x1x2 cells detailed earlier in this work. The lat- 
eral lattice constants are allowed to relax from a value of 
3.84 A at the bottom of the supercell (interfacial match- 
ing to Si(lOO)) to 3.99 A at the top of the supercell (the 
averaged unstrained lateral lattice constant of bulk or- 
thorhombic holmium silicide). The strain relief is just a 
linear function of height and, for example, we have not 
considered a case in which only the few layers near the 
interface are significantly strained. The strained island 
is shown from above in Fig. [7] 

Since we cannot currently fully optimise the c-axis in 
such a strained island we have simulated the blocking 



curves for four such islands, each with a different depth. 
This depth is varied by choosing two different c-axes for 
the two unit cells that form the island. These four c-axis 
values are; 1) both as bulk c-axis values (13.30 A), 2) both 
at the values obtained earlier in structure fitting with 
1x1x2 cells with no strain relief (13.14 A), 3) expansion 
of the c-axis of the cell nearest the interface by 2% and 
a bulk c-axis value in the top unit cell (13.57 A, 13.30 
A) and 4) expansion of the c-axis of the cell nearest the 
interface by 4% and the top ceU by 2% (13.83 A, 13.57 A). 
The calculated blocking curves are shown in Fig.[5]for the 
[lT0]/[110] geometry, the [111]/[211] geometry not being 
calculated because it would require averaging over two 
possible domains and thus twice as many simulations. 

On first viewing the incorporation of strain appears to 
make the fit worse (compare with the fit for the 1x1x2 
cell in Fig. 5(b)). The region with a scattering angle in 
the approximate range 117°-129° is especially bad. The 
scattering plane is shown schematically in Fig. [HI Wc can 
see from this that in the regions where the fit is bad it 
is the positions of the silicon atoms that are responsible 
for the poor fit of the blocking dips. Another minor sili- 
con atom blocking dip at around 107° is also a poor fit. 
Given that the positions and the thermal vibrations of 
the silicon atoms have not been optimised it is no sur- 
prise that the fit is made worse by the introduction of 
strain. 

However, the blocking dips produced by this structure 
are dominated by those caused by holmium atoms and 
close inspection of the major blocking dips at roughly 
82°, 113° and 135° supports the possibility of strain in 
these structures. The detail around these blocking dips 
is shown in Fig.[TOl In the strained cells as the c-axes are 
changed from the unrealistically low value obtained in 
the fit on an unstrained cell (13.14 A) through the bulk 
silicide value (13.30 A) and into the two cases where the 
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FIG. 8: Original in colour. Comparison of simulated strain in 
large 10x10x2 unit cells with experiment. The lateral lattice 
constant of each unit cell that makes up the nanoisland was 
allowed to relax from 3.84 A (interface lattice matching with 
Si(lOO)) to 3.99 A on the surface (bulk lattice constant). The 
four simulations plotted here correspond to fixed values of the 
two c-axes in the 10 x 10x2 unit cell. The first unit cell (cl) is 
nearest the surface and the second (c2) is nearest the interface. 



c-axis is expanded the blocking dip is shifted to become 
closer and closer to the experimentally measured dip. In- 
deed, at 84° and 113° the blocking dip for the strained 
cells is a better fit to experiment than that for the un- 
strained, fully optimised and multilayer averaged fit. 
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FIG. 9: Schematic slice showing the scattering plane for the 
[Il0]/[ll0] geometry and the atoms responsible for the princi- 
pal blocking dips observed. Filled circles represent Ho atoms 
and empty circles silicon atoms. 



Thus far, the question of the surface termination has 
not been addressed. The unit cells involved in the sim- 
ulations are holmium terminated which is not physically 
reasonable. Also, the LEED pattern shows a c(2 x 2) 
periodicity which is not accounted for in any of the simu- 
lated cells. We have attempted to investigate the surface 
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FIG. 10: Original in colour. Showing the detail around the 
major blocking dips in Fig. [S] With the gradual introduction 
of larger, more physically reasonable c-axis values into the 
strained cells the fit is improved until it is better than that 
for the fully optimised but non-strained cell. 



termination but the simulated blocking curves show very 
little sensitivity to an extra layer of silicon atoms that 
form a surface termination. Other experiments that are 
more sensitive to the very top atomic layer may be bet- 
ter suited to resolving the issue of the surface termination 
(such as MEIS itself at normal incidence with 50 keV He+ 
ions). 



V. LOWER COVERAGES 

The work presented here suggests that at a coverage 
of 6 ML and under our preparation conditions strained 
tetragonal silicide islands are formed when holmium is 
grown on Si(lOO). At lower coverages (<1 ML) nanowires 
form that are believed to be related to the hexagonal 
structure. The question naturally arises as to what struc- 
ture forms in the intermediate coverage regime. Does 
the tetragonal structure grow atop the hexagonal phase? 
Or is there a phase change at a certain minimum cov- 
erage which results in tetragonal rather than hexagonal 
growth? To attempt to answer these question we have 
taken experimental data at a coverage of 3ML under the 
same growth conditions. Fig. [TT] shows these data from 
the 3 ML sample compared with that from the sample 
grown using 6 ML of holmium. 

The principal blocking dips at 113° and 135° are 
present in both samples. In the 6 ML sample there are 
blocking dips at 90 , 95 and 100 that are not produced 
from the 3 ML sample. If we refer to Fig. [5] we can 
see that these blocking dips are produced by holmium 
atoms in the upper region of the structure. We would 
not expect to see these features in a thinner silicide layer. 
In the hexagonal structure the dominant blocking dip 
is at a scattering angle of approximately 95° (see fig. 
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FIG. 11: Original in colour. A comparison of the 6 ML cov- 
erage experimental data with data recorded after a sample 
was grown under identical growth conditions using a lower 
holmium coverage for the [ilO]/[110] geometry. The vertical 
scales have been artificially offset to aid the eye. 



4(a)). The simulations suggest that this feature should 
be present in hexagonal silicides that contain as few as 
3 layers of holmium. We can see no evidence for any 
blocking dips in this region in the experiment with 3ML 
coverage and we must conclude that even at this coverage 
the structure is tetragonal. 



VI. SUMMARY AND CONCLUSION 

Holmium silicide islands have been grown on the 
Si(lOO) surface and characterised using MEIS. Two struc- 
tures have been fitted to the experimental data and this 
process clearly shows that under our growth conditions 
it is the tetragonal phase of the silicide that is formed 
and not the hexagonal phase and that these two struc- 
tures do not coexist on the surface (the orthorhombic 
phase is too similar to the tetragonal phase for us to be 
able to determine if some of the orthorhombic nature of 
bulk HoSi is present in this surface tetragonal phase). A 
further experiment using a lower holmium coverage of 3 
ML has also been shown to have the tetragonal structure 
which confirms that the hexagonal phase is not the phase 
adopted at low coverages in this system. 

Other authors have reported growth of the hexagonal 
phase under similar growth conditions for some of the 
other RE silicides. It is interesting to speculate as to 
why similar experiments have reached different structural 



conclusions. It would appear that the structure that is 
formed is very sensitive to the growth conditions. Is- 
land morphology is very sensitively dependent upon the 
particular RE metal deposited and the annealing tem- 
perature used. In the paper by He et al. [22| they report 
a deposition rate of 0.5 ML per minute during the for- 
mation of DySi2 whereas Ye et al. [271 report 0.3 ML per 
minute. It could be that the mobility of the RE metal 
and/or the lattice mismatch anisotropy of the particular 
silicide in question are important enough factors during 
the very early stages of silicide formation to dictate the 
final structure. 

Using a simple 1x1x2 tetragonal unit cell the struc- 
ture fitting suggests a c-axis value that is too small to 
be physically reasonable when compared to the bulk 
structure. The blocking curves from a large 2500 atom 
10x10x2 nanoisland with lateral strain relaxation as a 
function of the distance from the interface have been 
simulated. The major blocking dips produced by this 
structure are a better fit to the experimental data than 
those produced by a non-strained periodic structure. The 
lesser blocking dips produced by silicon atoms in the sim- 
ulations do not show good agreement with experiment. 
This is to be expected because the thermal vibrations and 
the positions of these atoms have not been optimised in 
the structure fit. The results of the comparison of MEIS 
data with simulations clearly show that a physically rea- 
sonable c-axis value for the holmium silicide can only be 
obtained if the islands are strained to fit the Si(lOO) at 
the interface and the strain allowed to relax towards the 
top of the islands. 

A full structural fit for this system would require an op- 
timisation of the experiment-theory match with respect 
to all of the lateral and vertical spacings in the unit cell 
and their variation with distance from the interface, tak- 
ing into account lateral strain relief, vertical relaxation 
and the particular vibrations of each individual atom. 
This task is not currently tractable in terms of computa- 
tional resources and the available computer codes. 
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